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Abstract 

We present variational theory for optimal control over a finite time inter- 
val in quantum systems with relaxation. The corresponding Euler-Lagrange 
equations determining the optimal control field are derived. In our theory 
the optimal control field fulfills a high order differential equation, which we 
solve analytically for some limiting cases. We determine quantitatively how 
relaxation effects limit the control of the system. The theory is applied to 
open two level quantum systems. An approximate analytical solution for the 
level occupations in terms of the applied fields is presented. Different other 

applications are discussed. 
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The manipulation of quantum mechanical systems by using ultrashort time- dependent 
fields represents a challenging fundamental physical problem. In the last years, a consid- 
erable amount of experimental and theoretical work was concentrated on designing laser 
pulses having optimal amplitude and modulation. Thus the control of the quantum dynam- 
ics in various systems like atoms and molecules [I], quantum dots 0, semiconductors [[3], 
superconducting devices [§J and Bose-Einstein condensate f| was achieved. 

Several theoretical studies, most of them using numerical optimization techniques, have 
shown that it is possible to construct optimal external fields (e.g. laser pulses) to drive a 
certain physical quantity, like the population of a given state, to reach a desired value at a 
given time |] ||] . 

Although this kind of control might be relevant for many purposes, a more detailed 
manipulation of real systems may require the control of physical quantities over a finite 
time interval. The search for optimal fields able to perform such control is a much more 
challenging problem for which no theoretical description has been given so far. 

In this letter we present for the first time an analytical theory for the control of simple 
open systems over a finite time interval. By applying a variational approach we derive a 
high-order differential equation from which the optimal control fields are obtained. We also 
determine the influence of relaxation, the limits of this control and its potential applications 
to the manipulation of fundamental physical quantities, like the induced current through 
impurities in semiconductors or the population of electronic states at metallic surfaces. 

Our goal is to formulate a theory which permits to derive explicit equations to be satisfied 
by the optimal control field. Note that one can guess the form of such equations from general 
physical arguments. Since memory effects are expected to be important, one should search 
for a differential equation containing both the pulse area 9(t) = f to dt'V(t'), where V(t) is the 
external field envelope, and its time derivatives. Therefore, for the case of optimal control of 
dynamical quantities at a given time to, the differential equation satisfied by 9(t) must be of 
at least second order to fulfill the initial conditions 0(to), 0(t ). In the same way, the control 
of time averaged quantities over a finite time interval [to, to + T] with boundary conditions 
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requires a differential equation of at least forth order for 6{t) due to the boundary conditions 
for 8(t) and 9{t) at t and to + T. We show below that for certain open systems a forth order 
differential equation for the control fields arises naturally using variational approach as an 
Euler-Lagrange (EL) equation. 

We start by considering a quantum-mechanical system which is in contact with the 
environment and interacting with an external field E(t) = V(t) cos(ut). Here V(t) refers to 
an arbitrary pulse and u is the carrier frequency. The evolution of such system obeys the 
quantum Liouville equation for the density matrix p(t) with dissipative terms. The control 
of a time averaged dynamical quantity of the system requires the search for the optimal 
shape V(t) of the external field. 

Thus, in order to obtain the optimal V(t) on time interval [0, T] we propose the following 
Lagrangian (throughout the paper we use atomic units h=m=e=l) 

/T a r-T 

A(t)(-+iZ(t)}p(t)dt + (3 C x dt. (1) 

(3 is a Lagrange multiplier and A(t) is a Lagrange multiplier density. The first term in 
Eq. ([[]) ensures that the density matrix satisfies the quantum Liouville equation with the 
corresponding Liouville operator Z{t) |J. While the first term describes the dynamics of 
the system under the external field, the functional L\ explicitly includes the description of 
the optimal control and is given by 

dip, V) = C ob (p) + \V\t) + Ax (^r) \ (2) 

where A and Ai are Lagrange multipliers. £ b(p) refers to a physical quantity to be maximized 
during the control time. The second term represents a constraint on the total energy of the 
control field 

2 / E 2 (t)dt^ [ V 2 (t)dt = E . (3) 
Jo Jo 

The third term represents a further constraint on the properties of the pulse envelope. The 
requirement 
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dt < R, (4) 



where R is a positive constant, excludes infinitely narrow or sharp step-like solutions, which 
cannot be achieved experimentally. 

Assuming that the density matrix p(t) depends only on 9(t) and time, one obtains an 
explicit expression for the functional C\ = Ci(9,9,9,t). The corresponding extremum con- 
dition 5Ci = yields the high-order EL equation 

In order to solve Eq. (|5]) one can assume the natural boundary conditions 9(0) = 9(0) = 
6(T) = 0, 6(T) = &t, which also ensure that V(0) = V(T) = 0. The choice of the constant 
&t depends on the problem. In general, the constants 9t, R and Eq can be also object of 
the optimization. Note, that above formulated problem is highly nonlinear with respect to 
the function 9(t) and can be solved only numerically. 

Eq. (|5]) is the central result of this letter and provides an explicit differential equation 
for the control field. Note that this equation is only applicable if p = p(9(t),t). 

In order to show that Eq. @ can describe optimal control in real physical situations, we 
apply our theory to an open two level quantum system. This is characterized by the energy 
levels e± and 62, a dipole matrix element p and the longitudinal and transverse relaxation 
constants, 71 and 72, respectively. The carrier frequency of the control field is chosen to be 
the resonant frequency u — e% — e\. The dynamics of the density matrix pit) follows the 
equations (in the rotating wave approximation) 

.dp u 



dt 
.dpi2 



(-1) (pV(t)(p 2 i - P12) - 271P22), 

pV(t)(p 22 - Pu) ~ H2P12, (6) 



dt 

with £ = 1,2. Note that p n + p 2 2 = 1 and p 2 i — Pn- Eqs. (|6|) are used for the description 
of different effects, like for instance, the response of donor impurities in semiconductors to 
teraherz radiation 0, or the excitation of surface- into image charge states at noble metal 
surfaces [jlOfl . Therefore, the initial conditions are set as pxi = 1,P22 = Pn — P21 = 0. 



Eqs. (H) have the form idp(t)/dt = Z{t)p{t) and are difficult to integrate, since 
[Z(t), Z(t')] ^ 0. However, the commutators [Z(t), Z(t')] become arbitrarily small under 
the condition |Tf2 

d\ogV(t) 



-It 



« 1, (7) 



dt 

with £ = 1,2. In this case approximate solution for p 22 (t) is 

p 22 (t) = 2 ^(i)F- 1 (l - cosh(F) exp(-( 7l + T2 )t/2) 

+(7i + l2)t sinh(#) exp(-( 7l + 72 )t/2)#- 1 Y (8) 



where if = a/(( 7 i - 72 ) 2 £ 2 - 16 9 2 {t))/2, and F = 7i72 t 2 + 4 6> 2 (t). Note that this approx- 
imate solution becomes exact when 7l = 72 = or for a constant control field V(t) = Vq. 
The expression of Eq. (|j) has the form p = p(9(t),t) and therefore Eq. (|j) is applicable. 

Now we construct the functional C b(p) = P22W, so that the average occupation of 
the upper level n 2 = Jq P22(t)dt is maximized. Note, that n 2 proportional to the observed 
photocurrent M in teraherz experiments on semiconductors. The resonant tunneling current 



through an array of coupled quantum dots is also proportional to a such value [[□]] . 

We have calculated the optimal V(t) from the numerical integration of Eq. (|5|) for different 
values of the relaxation constants 7 i and 72 and of the energy Eq and the curvature R of 
the control fields. For simplicity we consider the control interval [0, 1]. 

In Fig. 1 we show the optimal field for an isolated ( 7 i = 72 = 0) and for an open two 
level system for given values of the pulse energy and curvature. Note, that for both cases 
the pulse maximum occurs near the beginning of the control interval. This leads to a rapid 
increase of the population p 22 (t) and therefore to a maximization of n 2 . In the case of an 
isolated system the pulse vanishes when the population inversion has been achieved, whereas 
for an open system the pulse must compensate the decay of p 22 (t) due to relaxation effects 
and remains finite over the whole control interval. 

In the inset of Fig. 1 we show the corresponding dynamics of the population p 22 (t) for 
both cases. As mentioned before, Eq. (|8|) is exact for the isolated system. Note, that for the 



open system the analytical form of p 22 (Eq. @) compares well with the numerical solution 
of the Liouville equation. This indicates that V(t) fulfills the condition ([7]) on the control 
interval. 

We found that the value of n 2 increases both for the isolated and for the open system 
monotonously with energy of the optimal control field. In Fig. 2 we plot n 2 as a function 
of the energy Eq and the curvature R of the optimal fields obtained from Eq. (f|). Note, 
that pulses of fixed shape (for instance Gaussian) would show an oscillating behavior for 
increasing energy due to Rabi oscillations [f|. The monotonous increase is a feature which 
characterizes the optimal pulses. 

In order to achieve a simplified study of the physics contained in the control fields of 
Fig. 1, we analyze the problem in certain limiting cases. For instance, if jipT <C 1 one 
can neglect decoherence within the control interval and Eq. (^) becomes p 22 (t) — sin 2 (8(t)). 
In order to make the problem analytically solvable, we reduce the order of the differential 
equation for the control fields. For that purpose we replace the constraint on the derivative 
of the field envelop (Eq. (^)) by a weaker one obtained from the condition 



jf e 2 (t)dt-^J^ e(t)dt) = J c w (0)>s, 



(9) 



where S is a positive constant. Eq. @ merely bounds the width of the envelope V(t) in 
order to avoid unphysically narrow pulses. Thus, the Lagrangian density Ci for the optimal 
control has the form 

C 1 =p 22 (t) + \6 2 (t) + \ 2 C w (6), (10) 

while the corresponding EL equation is given by 

2\'§(t) - sin(2 0(t)) = 0. (11) 

Note, that condition (||) only leads to a rescaling of the Lagrangian multiplier A to A' = A+A2. 
The second order differential Eq. (|TT|) requires two boundary conditions, for which we choose 
6(0) = and 9{T) = ir/2 (which ensure the population inversion). Eq. ( |TT]) resembles that 
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for a mathematical pendulum and can be solved analytically. The resulting field envelope 
is given by 

V(t) = V(p)dn(V(0)t,C), (12) 

where dn is the Jacobian elliptic function, and C = — (\'V 2 (0)) 1 is a constant of integration. 
Note, that V(0) ^ 0. Using conditions (f|) and @ we determine coefficients A and A2. If 
we choose C — > 1 then we can obtain V(T) — > 0. In this case, Eq. fll2|) can be significantly 
simplified to V{t) = f arccos[2exp (V{0)t)/{1 + exp (2V{0)t))]. 

In Fig. 3 we plot the optimal control field V(t) which maximizes the Lagrangian ([T0| ) for 



isolated and open two level systems. In both cases the field has its maximum value at t — 
and exhibits a monotonous decay. As in the case of the solutions of the forth-order Eq. ([5]) 
the control field is broader for the open system. In the inset of Fig. 3 we plot the population 
P22(t). The overall behavior of is similar to that of the populations shown in Fig. 1. 

It is important to point out, that a Lagrangian of the form of Eq. ( |T0D always leads to a 
second order differential equation for the control fields as long as the condition p = p(8(t),t) 
is satisfied. Therefore, one cannot demand extra boundary conditions for the fields V(0) = 
V(T) = 0. Otherwise one would obtain the trivial solution V(t) = 0, which is not consistent 
with either (||]) or (||]). Therefore, if conditions on V(0) and V(T) have to be imposed, a 
Lagrangian leading to a forth order differential equation is necessary, as we have shown 
before. 

As it was mentioned before increases monotonously with the pulse energy for the 
optimal field. Since for the isolated system n 2 approaches the maximum possible value 
i%2 = T, in the case of nonisolated systems there is a limit. In order to show that this limits 
is due to general physical reasons we analyze the occupation p22(t) (Eq. (|8|)) in more detail. 
For a strong control field satisfying ji^t/Oft) 1 the occupation P22W always lies under 
the curve p™ 2 ax {t) — (1 + exp( — (71 + 72)t/2))/2. This means that it exhibits an absolute 
upper bound. Therefore due to dissipative processes the following inequality holds for the 
controlled averaged value of P22: 



na= / P22(t)^<T/2 + (l-exp(-(7 1 + 7 2 )T/2)/( 7l + 7 2 ). (13) 

Eq. fll3l) shows the absolute limit for the optimal control of averaged occupations in 
open two level systems. In Fig. 4 we show the maximal possible value p™^ x (t) and the time 
evolution of P22W induced by 40 randomly generated pulses (for some of which the condition 
li,%t/Q(t) *C 1 is even not strictly fulfilled). From Fig. 4 we conclude that under the action 
of arbitrary control fields, the life time of the upper level cannot be longer than 2/ (71 +72). 

Using this result we can determine the maximal possible life-time for an image state at 



a Cu(lll) surface which can be achieved by pulse shaping. According to Hertel et al. [10 
those states are characterized by 71 = 5 ■ 10 13 s _1 and 72 = 7i/2. Thus, our theory predicts 
in that case an effective decay constant j e ff = (71 + 72V2 = 3.75 • 10 13 s _1 . 

In summary, we presented a theory for the description of optimal control of time-averaged 
quantities in open quantum systems. In particular we have shown that the boundary con- 
ditions of the problem make a significant influence on the shape of the optimal fields. In 
contrast to other approaches our theory allows to derive an explicit differential equation for 
the optimal control field, which we integrated both numerically and exactly for some limiting 
cases. Our approximation p(t) = p{6, t) was checked by direct integration of the Liouville 
equations and it seems to hold also in the case of strong relaxation. Using our theory we 
found the optimal fields which maximize the population of the upper levels of isolated and 
open two-level systems. We found an absolute upper bound for this kind of optimal control. 
Our approach can be used for further investigations, for instance, control of the dynamics of 
multi-level systems. 
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Time (a.u.) 

FIG. 1. Optimal control field for an isolated two level system (71,2 = 0, solid line). The pulse 
energy is Eq = 4.57 and the pulse curvature R = 128.4. The dashed line shows the optimal pulse 
for the open system (71 = 272 = 5) with energy Eq = 53.54 and curvature R = 808.8. Inset: 
Dynamics of the occupation P22(t) for an isolated system (thick solid line) and with relaxation 
(dash dotted line- using formula (H), thin solid line-numerical solution of the Liouville equation 
@). Arbitrary units are used. 
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FIG. 2. Dependence of the averaged occupation n2 as a function of the energy Eq and the 
curvature it! (see Eq. (||) ) of the optimal pulses. 
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FIG. 3. Optimal control field for an isolated two level system (71,2 = 0, solid line). The pulse 
energy is Eq = 20.50. Dashed line: optimal field for the open system (71 = 272 = 5) with a pulse 
energy Eq = 89.72. Inset: Dynamics of the occupation p22(t) f° r an isolated system (thick solid 
line) and with relaxation (dash dotted line using Eq. (g), thin solid line- numerical solution of the 
Liouville equation (ph). 
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FIG. 4. Dynamics of the occupation P22(t) for 40 randomly generated control pulses and for 
7i = 272 = 1 (thin solid lines). The thick solid line represents a bound for the possible values of 
ffi*(t) = (1 + exp(-( 7l + 72 )t/2))/2. 
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